Hydrodynamic Pair Attractions Between Driven Colloidal Particles 
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Colloidal spheres driven through water along a circular path by an optical ring trap display 
unexpected dynamical correlations. We use Stokesian Dynamics simulations and a simple analytical 
model to demonstrate that the path's curvature breaks the symmetry of the two-body hydrodynamic 
interaction, resulting in particle pairing. The influence of this effective nonequilibrium attraction 
diminishes as either the temperature or the stiffness of the radial confinement increases. We find 
a well defined set of dynamically paired states whose stability relies on hydrodynamic coupling in 
curving trajectories. 
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Position correlations among particles in thermal equi- 
librium stem solely from their potential energy. Although 
the motions of fluid-suspended particles are correlated 
via flows that they induce (hydrodynamic interactions), 
this has no effect on the equilibrium configurations. For 
particles driven out of equilibrium, however, this is no 
longer the case, and they may get closer or drift apart 
in the absence of an interaction potential. Examples are 
found in studies of nonequilibrium particle aggregation 
and droplet coalescence under flow (e.g., (H), aggregation 
of blood cells in channel flow|2[ , and strong correlations 
in sedimcnting suspensions [3(. Nonequilibrium pair at- 
tractions have been reported in simulations of particles, 
embedded in a flowing bath of smaller particles or 
driven through such a bath Yet, all these phenomena 
are observed in the presence of cither many-body effects 
or nonuniform flow, where spatially separated particles 
experience a different drive. 

It is not self-evident that two identical particles, having 
no intrinsic interactions, can attract or repel each other 
under conditions of identical driving jf|. This is because 
a uniaxial driving force, breaking rotational symmetry, 
is generally insufficient to break the symmetry of the hy- 
drodynamic interaction on the two-body level. A famil- 
iar example is an isolated pair of noninteracting colloidal 
particles, sedimcnting under gravity in an unbounded vis- 
cous fluid. One particle drags the other to the exact same 
extent that the other pushes it; hence, their separation 
remains constant. One way to break this symmetry is 
to introduce a boundary. Thus, those two sedimcnting 
particles, when falling away from a horizontal wall, ex- 
erience a hydrodynamic force that draws them together 
Another way is to consider more than two particles. 
For example, driving a one-dimensional array of particles 
by a uniform flow through a slit creates hydrodynamic 
restoring forces, resulting in density waves [7[. 

In this Letter we show that the symmetry of the hydro- 
dynamic interactions of identically driven particles in a 
fluid can be spontaneously broken on the two-body level 
if their trajectories curve. The result is the formation 



of dynamically paired states, whose stability determines 
the steady transport properties of the system. 

Our samples consist of aqueous dispersions of colloidal 
polystyrene particles, sealed in a 50 /xm thick gap be- 
tween a glass microscope slide and a coverslip. The 
Debye-Hiickcl screening length is of order 10 nm, imply- 
ing that the charge-stabilized spheres' electrostatic in- 
teractions are negligible at micrometer-scale separations. 
Individual spheres are captured and driven through water 
by extended optical traps, created with the holographic 
optical trapping technique [H, In our implementa- 
tion [T3|, we use a helical mode of light, which focuses 
to a ring-like optical trap known as an optical vortex 
@, 3 HI ■ Particles arc drawn to the optical vortex by 
forces directed by intensity gradients [12| |. and are pro- 
pelled around it by radiation pres sure arising from the he- 
lical beam's phase gradients [13)] • Adaptive optics meth- 
ods are employed to reduce variations in the tangential 
driving force [9|, [l4|. The optical force profile, as inferred 



from the motion of a single particle, is shown in Fig.[TJa), 
revealing a roughness of about 30%. The optical vortex is 
powered with 2.5 W of laser light at a vacuum wavelength 
of 532 nm, which suffices to drive trapped particles with 
tangential velocities Vg of a few tens of micrometers per 
second, corresponding to angular velocities Q of a few 
radians per second. Standard methods of digital video 
microscopy [l5T ] are used to track the particles at a rate 
of 30 frames/s with 20 nm accuracy. 

We present results for two sets of experimental param- 
eters (particle diameter a, vortex radius R, tangential 
driving force Fg, and radial trapping stiffness fc r ), as 
listed in Table HI Experiments were performed for num- 
bers of particles ranging from N = 1 to N = 12, running 
around the ring in single-file trajectories. The particles' 
inertia may be ignored because of the negligible Reynolds 
number, Re = pVgcr/ri ~ 10~ 5 , where r\ and p are the vis- 
cosity and mass density of water, respectively. Random 
thermal forces contribute fluctuations to the circulating 
particles' velocities [l6| . The strength of the driving rela- 
tive to thermal fluctuations is characterized by the Peclet 
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FIG. 1: (color online) (a) Magnitude of the tangential com- 
ponent of the force exerted by an optical vortex (system A) 
on a single colloidal sphere, (b) Visualization of a pair having 
an angular separation AO. (c,d) Distributions of experimen- 
tally measured angular separation (c) and ratio between the 
radii of the leading and trailing particles (d) in a two-particle 
vortex for the two systems. 



number, Pc = Fga / (k-gT) ~ 90 and 74, respectively, for 
systems A and B, where fee is Boltzmann's constant, and 
T the absolute temperature. The ratio of the radial con- 
finement strength and tangential driving, n = k r R/Fg, 
serves as another characteristic parameter. We estimate 
k — 30 for both systems. 

Because the particles have no intrinsic long-range in- 
teractions and because they are driven by the same op- 
tical force, a pair of particles might be expected to dif- 
fuse independently and thus sample all possible separa- 
tions with uniform probability. Experimentally, however, 
the two particles enter into a long-lived paired state [l7| , 
as quantified in the histograms of angular separations, 
Fig- He)- The narrower distribution observed for system 
A is consistent with its larger Pe. 

The observed correlation between the particles' mo- 



system 


a (/xm) 


R (fjxo.) 


Fe ( P N) 


k r (pN//im) 


A 


1.48 ± 0.08 


6.25 ±0.05 


0.25 ±0.08 


1.2 ±0.3 


B 


0.97 ±0.03 


6.19 ±0.05 


0.31 ±0.09 


1.7 ±0.8 



TABLE I: The two experimental parameter sets: particle di- 
ameter a, vortex radius R, tangential driving force F$, radial 
confinement stiffness k r . We calculate Fe from the mean tan- 
gential velocity in a single-particle vortex, Fg = Ve/B s , and 
k r from the single particle's mean-square radial fluctuations 
assuming thermal equipartition [l(J, k r — ksT / ((At) 2 ) . 



tions must arise from the interplay of optical and hydro- 
dynamic forces acting on them. A single sphere circulates 
around a smooth optical vortex with a uniform tangential 
velocity, Vg = B s Fg, where B s = (3nr]a) _1 is the sphere's 
Stokes self-mobility. As a sphere moves, it creates a flow 
field G(r)Fg at position r. A second sphere placed on the 
optical vortex at r experiences not only the optical force 
Fg, but also the drag force due to its neighbor's flow. 
To leading order in <r/r, the change in the tangential 
velocity of the second particle due to its neighbor is sim- 
ply the local tangential flow velocity, and G is given by 
the Oseen tensor [l8[, Gij(r) = (87r77r) _1 (6ij + rjrj/r 2 ), 
i,j = x,y,z. Applying this result to our case of two 
spheres on a ring of radius R, separated by an angle 
9ij = 6j — 0i , we get for the tangential velocity, 
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Equation (TTJ is symmetric under particle exchange, 9ij — > 
8ji = —6ij, implying that the tangential component of 
the hydrodynamic coupling cannot account for the ob- 
served pairing. Because the circulating particles follow 
a curved path, however, the flow field due to particle i 
generally has a component normal to the trajectory at 
the position of particle j. (See [l?} for a theoretical vi- 
sualization of the flow field.) Since the particle is held 
radially by the trap, this normal flow component gives 
rise to a radial force, 
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This hydrodynamic force displaces the particle radially 
until it is counter-balanced by an equal optical force. 

The fact that h(9ij) is antisymmetric under particle 
exchange implies that the two particles are displaced op- 
positely in the radial direction. As is depicted schemati- 
cally in Fig. [21 the leading particle is shifted to a slightly 
larger radius, Ri > R, and the trailing particle to a 
slightly smaller one, f?2 < R- This spontaneous sym- 
metry breaking is observed in the statistics of measured 
trajectories, plotted in Fig. [ljd). If the tangential force 
Fg is assumed constant, the radial displacement increases 
the trailing particle's angular velocity relative to that of 
the leading particle, and so causes the particles to form 
a hydrodynamically bound pair. 

Pairing also appears in the dynamics of systems com- 
prising more than two spheres. Driven rings of hydrody- 
namically coupled spheres thus constitute an interesting 
model system for studying spontaneous symmetry break- 
ing in many-body systems out of equilibrium. FigureEta) 
shows the measured collective mobility as a function of 
particle number N. The increase in the circulation rate 
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FIG. 2: (color online), (a-c) The pairing mechanism. As 
optical forces drive the particles along the ring, the resulting 
fluid flow pushes the leading particle outward and pulls the 
trailing one inward. Once displaced radially, the particles' 
angular velocities change, making the trailing particle catch 
up with the leading one. (d,e) Simulation results for the effect 
of radial confinement strength (d) and temperature (e) on the 
average angular separation (green circles) and average angular 
velocity (red squares) in a two-particle vortex. The insets 
show the distributions of radial positions of the leading (Ri) 
and trailing (-R2) particles, and the pair's angular separation, 
for T = 0.0375 and k = 25. 



due to a decrease in the hydro-dynamic drag is more pro- 
nounced for even values of TV. In systems with odd TV 
the tendency to form pairs leaves one particle unpaired, 
which reduces the influence on the collective drag. 

Let us focus on the case of even N and consider an ide- 
alized configuration of N/2 pairs, evenly positioned along 
the ring. In the absence of hydrodynamic interactions, 
each pair will move with tangential velocity Vg = B p Fg, 
where B p = aB s is the mobility of a pair (a being a 
prefactor of order 1). Including tangential hydrodynamic 
couplings as above, we get for the collective mobility, 
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where g(9) has been defined in Eq. (jTJ) . In Fig. Eta) we 
compare the experimentally measured collective mobil- 
ity, for even N, with the one predicted by Eq. ((31), using 
the vertical shift a as a fitting parameter. Despite the 
simplifying assumptions underlying Eq. ([3]), a qualitative 
correspondence is obtained for a = 1.10 ± 0.02. We at- 
tribute the larger discrepancies for N = 2 and N = 4 to 
larger deviations from the idealized configuration of fixed 
pairs in these dilute systems. 

We systematically investigated how the pairing dy- 
namics depend on Pe and K by performing Stokesian 
Dynamics simulations [III]. The simulated particles are 
restricted to the plane and confined to a ring of radius 
R = 5cr by a harmonic radial potential of stiffness k r . A 
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FIG. 3: (color online). Collective mobility, scaled by the 
single-particle mobility, as a function of particle number, (a) 
Experimental results for even N (green circles) and odd N 
(blue squares), measured from the ratio of mean velocities for 
N particles and for a single particle in system A. (b) Simula- 
tion results for even N (green circles, T = 0.1; green triangles 
T = 0.02) and odd N (blue squares, T = 0.1; blue diamonds, 
T — 0.02), at k = 25. The solid lines are theoretical fits 
for even N [Eq. ©] with a = 1.10 ± 0.02 (experiment) and 
a = 1.57 ± 0.05 (simulation). Error bars are smaller than 
symbols unless otherwise indicated. 



constant tangential force Fg drives the particles around 
the ring, and a repulsive Weeks-Chandlcr-Anderson pair 
interaction [i(| accounts for collisions. Hydrodynamic in- 
teractions between particles are modeled using the many- 
body Rotne-Prager mobility tensor [2l|,[22j]. Temperature 
is introduced via a Gaussian-distributed random force 
that obeys the fluctuation-dissipation relation using the 
same mobility tensor. The angular velocity is measured 
in units of B s Fg/R and the simulation temperature is 
measured in units of Fgo 1 /(fes-R). In these units, the 
Peclet number is Pc = (Rfa)T- 1 = 5/T. 

Figure [UJd) presents simulation results for two parti- 
cles on a ring, confirming the role of radial symmetry 
breaking in creating particle pairs. At low temperature 
and weak radial confinement, the particles form perma- 
nent pairs. Paired particles follow trajectories of differ- 
ent radii with R1/R2 ~ 1.01 [Fig. [U[d) inset]. Pairing in 
simulation is more pronounced than in experiment, with 
a smaller and more narrowly distributed angular sepa- 
ration [Fig. (2Ie) inset] [23}. For infinitely strong radial 
confinement, no pairing is observed. Vigorous diffusion 
at elevated temperatures overcomes the hydrodynamic 
pairing mechanism and allows pairs to move apart. This 
increases their drag, and reduces their mean angular ve- 
locity [Fig. He)]. 

Simulations for rings with even N > 2 exhibit pairing 
as well. Comparing the collective mobilities with our sim- 
plified model [Fig. H^b)], we obtain a reasonable fit using 
the value a = 1.57 ± 0.05. The value of a = B p /B s can 
be related also to the inter-particle distance through the 
theoretical expression for the mobility of a pair of spheres 
[Til ] . In the simulation we measure an angular separation 
of A8 = 12.9 ± 0.1°, which yields a = 1.52 ± 0.02. These 
two independent measurements of the pair mobility agree 
to within the statistical error. The smaller experimental 
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value of a [Fig. EJa)] suggests still weaker pairing. These 
quantitative differences between experiment and simula- 
tion may be related to such effects as the roughness of 
the experimental force landscape [l6j and a possible non- 
intrinsic repulsion (23j . 

At low temperatures, particles may assemble into dif- 
ferent configurations as a result of their tendency to form 
pairs, giving rise to a rich behavior of degenerate limit 
cycles. The simplest example is the case of three parti- 
cles The breaking of a triplet into a faster pair 
and a slower single particle exists also in the absence of 
curvature and radial freedom, yet, in simulations we ob- 
serve long-term pairing and the consequent limit cycle of 
an alternating pair only when the radial confinement is 
made sufficiently weak. 

A simulated four-particle vortex is found to converge 
into three stable limit cycles, shown in Fig. |4j For a 
given radial confinement strength, particles adopt con- 
figurations 1 and 3 at low temperatures and configura- 
tions 2 and 3 at higher temperatures. With increasing 
temperature the pairs gradually break. In configuration 

1 the two pairs are separated by a certain angle (0.67T in 
our simulation), which seems unrelated to the system's 
geometry. The transition between configurations 1 and 

2 explains the difference in collective mobility seen in 
Fig. EJb). Since configuration 2 is compatible with our 
calculation of equidistant pairs, we get a better agree- 
ment for the higher temperature in which it is stabilized. 
This effect might also contribute to the increased mobil- 
ity of four particles seen experimentally [Fig. EJa)]. The 
assumed experimental Peclet numbers [16J correspond to 
simulation temperatures of T = 0.055 and 0.068 for sys- 
tems A and B, respectively, with n values around 30. 
This might indicate that the experiments lie in the inter- 
mediate region between pairing and unpairing (Fig. |4]) . 
As the number of particles increases so does the number 
of limit cycles, leading to a more complex picture to be 
presented elsewhere. 
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FIG. 4: Simulation results for the normalized average an- 
gular velocity Q and particle configuration, as a function of 
temperature, in a four-particle vortex. Insets illustrate the 
three configurations. 

The spontaneous symmetry breaking and dynamic 



coupling presented here can arise whenever particles are 
driven through a fluid medium along curving trajectories, 
provided that they are allowed to slightly move in the 
radial direction. Such nonequilibrium pair attractions, 
therefore, should affect the behavior of diverse out-of- 
equilibrium colloidal systems. In addition, we intend to 
use this model system to study further issues of nonequi- 
librium physics in a well controlled environment. 
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